# importação das bibliotecas Numpy, Matplotlib e Seaborn e do interpretador interativo.
import warnings
warnings.simplefilter(action = 'ignore', category = FutureWarning)
import numpy as np
import seaborn as sns
import matplotlib.pyplot as plt 
import matplotlib.ticker as mticker
sns.set(style = 'darkgrid')
from matplotlib.animation import FuncAnimation
from IPython.display import HTML

# Vamos escolher três valores diferentes de volume de gás para calcular as trajetórias dos processos.
# O ponto 1 = (P1, V1) define o início do ciclo. Os valores e as unidades adotados são arbitrários e
# faremos nR = 1, por simplicidade.
V1 = 2.0
V2 = 1.0
V3 = 2.0
V4 = 4.0

# Ponto 1: Início do primeiro processo - compressão isotérmica do gás. 
# No Ponto 1 do ciclo, o volume do gás é V1 e a pressão é P1:
# Podemos encontrar o valor da pressão P2 usando os valores dos volumes V1 e V2, conhecendo P1.
T1 = 2.0  # Temperatura no Ponto 1 do Diagrama PV.
P1 = T1 / V1  # Pressão no Ponto 1 do Diagrama PV.
# Ponto 2: Final do primeiro processo - compressão isotérmica do gás.
# Gerando os pontos para a curva que representa a compressão isotérmica do gás.
V_comp_isotermica = np.linspace(V1, V2, 100)
P_comp_isotermica = P1 * (V1 / V_comp_isotermica)

# Ponto 2: Final do primeiro processo - expansão isotérmica do gás.
T2 = T1 
P2 = P1 * (V1 / V2)
# Ponto 2: Início do segundo processo - expansão isobárica do gás. 
P3 = P2
# Gerando os pontos para a curva que representa a expansão isobárica do gás.
V_exp_isobarica = np.linspace(V2, V3, 100)
P_exp_isobarica = np.linspace(P2, P3, 100)

# Ponto 3, final do segundo processo - expansão isobárica do gás. Vamos usar a equação 
# de estado dos gases ideais: PV = nRT, com nR = 1, por simplicidade. Logo: PV = T.
# No Ponto 3 do ciclo, o volume do gás é V3 e a pressão é P3.
T3 = T2 * (V3 / V2)
# Ponto 3, início do terceiro processo - expansão isotérmica do gás. 
P4 = P3 * (V3 / V4)
# A temperatura no Ponto 4 é igual à temperatura no Ponto 3 (isoterma):
T4 = T3 
# Gerando os pontos para a curva que representa a expansão adiabática do gás.
V_exp_isotermica = np.linspace(V3, V4, 100)
P_exp_isotermica = P3 * (V3 / V_exp_isotermica)

# Ponto 4, final do terceiro processo - expansão isotérmica do gás.
# Ponto 4, início do quarto processo - compressão isobárica do gás.
P1 = P4
# Gerando os pontos para a curva que representa a compressão isobárica do gás.
V_comp_isobarica = np.linspace(V4, V1, 100)
P_comp_isobarica = np.linspace(P4, P1, 100)

# Plotagem do diagrama PV.
plt.subplots(1, 1, figsize = (10, 6), dpi = 200)
plt.plot(V_comp_isotermica, P_comp_isotermica, 'orange', label = 'Compressão Isotérmica')
plt.plot(V_exp_isobarica, P_exp_isobarica, 'b', label = 'Expansão Isobárica')
plt.plot(V_exp_isotermica, P_exp_isotermica, 'r', label = 'Expansão Isotérmica')
plt.plot(V_comp_isobarica, P_comp_isobarica, 'g', label = 'Compressão Isobárica')
plt.xlabel('Volume (U.A.)')
plt.ylabel('Pressão (U.A.)')
plt.xlim((0, 5.0))
plt.ylim((0.5, 2.5))
plt.title('Ciclo Ericsson no Diagrama PV')
plt.legend()
plt.grid(True)
# Setas indicando os inícios dos processos do ciclo.
plt.annotate('1', (V1, P1), (V1 + 0.32, P1 + 0.32), arrowprops = dict(arrowstyle = '->', color = 'black'), fontsize = 12)
plt.annotate('2', (V2, P2), (V2 + 0.3, P2 + 0.3), arrowprops = dict(arrowstyle = '->', color = 'black'), fontsize = 12)
plt.annotate('3', (V3, P3), (V3 + 0.3, P3 + 0.3), arrowprops = dict(arrowstyle = '->', color = 'black'), fontsize = 12)
plt.annotate('4', (V4, P4), (V4 + 0.32, P4 + 0.32), arrowprops = dict(arrowstyle = '->', color = 'black'), fontsize = 12)
plt.gca().yaxis.set_major_locator(mticker.MultipleLocator(0.5))
plt.show()

import matplotlib.pyplot as plt
from matplotlib.animation import FuncAnimation
from IPython.display import HTML

plt.rcParams['animation.embed_limit'] = 256
fig, ax = plt.subplots(figsize=(10, 6))
ax.set_ylim(0, 5)
ax.set_xlim(0, 5) 
ax.set_xlabel('Volume (U.A.)')
ax.set_ylabel('Pressão (U.A.)')
ax.set_title('Animação Ericsson no Diagrama PV')

V_A = V_comp_isotermica
P_A = P_comp_isotermica
V_B = V_exp_isobarica
P_B = P_exp_isobarica
V_C = V_exp_isotermica
P_C = P_exp_isotermica
V_D = V_comp_isobarica
P_D = P_comp_isobarica



# Inicializando as curvas como linhas vazias
line_A, = ax.plot([], [], 'orange', label = 'Compressão Isotérmica')
line_B, = ax.plot([], [], 'b', label = 'Expansão Isobárica')
line_C, = ax.plot([], [],  'r', label = 'Expansão Isotérmica')
line_D, = ax.plot([], [], 'g', label = 'Compressão Isobárica')


# Adicionando a legenda e a grade
ax.legend()
ax.grid(True)

# Função de inicialização, chamada apenas uma vez
def init():
    line_A.set_data([], [])
    line_B.set_data([], [])
    line_C.set_data([], [])
    line_D.set_data([], [])
    return line_A, line_B, line_C, line_D

# Função de animação
def animate(frame):
    # Atualizando os dados das curvas
    if frame < len(V_A):
        line_A.set_data(V_A[:frame], P_A[:frame])
    elif frame < len(V_A) + len(V_B):
        line_B.set_data(V_B[:frame-len(V_A)], P_B[:frame-len(V_A)])
    elif frame < len(V_A) + len(V_B) + len(V_C):
        line_C.set_data(V_C[:frame-(len(V_A)+len(V_B))], P_C[:frame-(len(V_A)+len(V_B))])
    else:
        line_D.set_data(V_D[:frame-(len(V_A)+len(V_B)+len(V_C))], P_D[:frame-(len(V_A)+len(V_B)+len(V_C))])

    return line_A, line_B, line_C, line_D

# Configurando a animação
total_frames = len(V_A) + len(V_B) + len(V_C) + len(V_D)
ani = FuncAnimation(fig, animate, frames=total_frames, init_func=init, blit=True, interval=100)

# Exibindo a animação usando JSAnimation
HTML(ani.to_jshtml())

Ciclo de Ericsson¶